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Dark Matter (DM) may be a thermal relic that annihilates into heavier states in the early Universe. 
This Forbidden DM framework accommodates a wide range of DM masses from keV to weak scales. 
An exponential hierarchy between the DM mass and the weak scale follows from the exponential 
suppression of the thermally averaged cross section. Stringent constraints from the cosmic microwave 
background are evaded because annihilations turn off at late times. We provide an example where 
DM annihilates into dark photons, which is testable through large DM self-interactions and direct 
detection. 


Introduction: Dark Matter (DM) accounts for over 
80% of the matter of our Universe, but its particle ori¬ 
gin remains a mystery. If it is a thermal relic, its energy 
density today depends on its annihilation rate, (cru), in 
the early Universe [T], 
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The observed DM abundance, ^ 0.1, is repro¬ 

duced if DM annihilates with a weak scale cross section. 

Typically, (av) ^ where is the DM mass 

and ad characterizes its interaction strength. For siz¬ 
able ad, the DM mass should be near the weak scale, 
m,/, ~ 0.1 — 10 TeV. Direct detection experiments have 
been making rapid progress searching for weak scale DM 
scattering against nuclei, and significant reach is pro¬ 
jected over the next decade (see Ref. [2] for a review). 
However, pending a possible discovery, it is important to 
cast a wide net by exploring a diverse set of models with 
varying DM masses and signatures. 

There are several classes of models where the DM mass 
may be much lighter than the weak scale. One possibility 
is to take ad ^ 1, keeping c?Jrv?^ fixed [2- A second 
possibility is that the DM energy density is determined 
by an asymmetry, similarly to baryons, such that the DM 
mass varies inversely with the size of the asymmetry [U 
[ 3 . A third possibility is that the DM relic density is 
determined by 3 —>■ 2 annihilations instead of 2 —?> 2 
annihilations, requiring light DM to achieve the necessary 
cross section [SE!. 

In this letter we propose a novel framework with nat¬ 
urally light DM. We assume that DM dominantly anni¬ 
hilates into heavier particles: iliip —> xy with A- my > 
2m,/,. Such annihilations are called forbidden channels 
because they vanish at zero temperature. Forbidden 
channels can proceed at finite temperature in the early 
universe, due to the thermal tail with high velocity 
The thermally averaged cross section for forbidden chan¬ 
nels is exponentially suppressed, as we derive below. This 
exponential suppression allows for a weak scale cross sec¬ 
tion when a^/m^ is much larger than the weak scale, 
permitting light m.^ with a sizable ad- While it is an 


old observation that forbidden channels can dominate 
DM annihilations [S] , previous work has focused on weak 
scale DM (see for example Refs. [SHU]) and metastable 
relics m- Forbidden channels remain unexplored as a 
mechanism for light DM. 

We call Forbidden Dark Matter the class of models 
where the DM relic abundance is dominantly set by for¬ 
bidden channels. Models of light DM are highly con¬ 
strained by the Cosmic Microwave Background (CMB), 
which is sensitive to energy injection into the photon 
plasma at the recombination epoch [Ml [15]. Forbid¬ 
den DM evades these constraints because annihilations 
shut off at low temperatures. Tree ^ mx + my — 2m,/,. 
This framework naturally includes large self-interactions, 
ipip —> ifij}, which are unsuppressed and therefore have 
an exponentially larger cross section than the forbid¬ 
den annihilation rate. Large DM self-interactions may 
help address problems with structure formation in non¬ 
interacting DM scenarios [HM] and may be indicated 
by recent cluster observations [52 . 

For the remainder of this letter, we consider an ex¬ 
ample model that illustrates the key features of the For¬ 
bidden DM framework. We take the DM to be a Dirac 
fermion, if, that is neutral under Standard Model quan¬ 
tum numbers and charged under a hidden U{l)d gauge 
symmetry that is broken at low energies. DM annihilates 
into hidden photons, tpip —)■ ')dld (see the hrst diagram in 
Fig. [^. The parameters of the model are (m,/,, m-^^, a^). 
We consider the forbidden portion of parameter space, 
m,/, < m^,^. The non-forbidden parameter space of this 
model was studied by Refs. [2^ 125] . 

The organization of the rest of this letter is as follows. 
First, we discuss the computation of the relic density. 
Second, we derive the DM self-interaction cross section. 
Third, we consider the possibility that the hidden pho¬ 
ton couples to the Standard Model (SM) through gauge 
kinetic mixing, leading to signals in direct and indirect 
detection. Finally, we consider the possibility that the 
hidden sector is thermally decoupled from the SM, al¬ 
lowing for sub-MeV DM masses. We finish with our con¬ 
clusions. In an upcoming paper, we will consider light 
DM from co-annihilations, where an exponential factor 
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FIG. 1. The left panel contains Feynman diagrams relevant for (1) the relic density, (2) self-interactions, (3) indirect detection, 
and (4) direct detection. The right panel shows the relic density, h?, as a function of the mass splitting A = (m.yj — . 

The red (blue) curves correspond to = 1 GeV (MeV) and the solid (dashed) curves correspond to cxd = 0.1 (10“®). 


also allows for very light DM |26) . 

Relic Density: The relic density of Forbidden DM is 
determined by the solution of its Boltzmann equation, 

+ 3Hn^ = - v) v) , ( 2 ) 


where &re the number densities, H is the Hubble 

parameter, v) denotes the thermally averaged 

(inverse-)annihilations, and we have assumed that 7 ^ re¬ 
mains in equilibrium during freeze-out. The solution is 
approximately given by Eq. ([^, with the annihilation 
rate given by For simplicity, Eq. (j^ neglects 

the dependence on the number of relativistic degrees of 
freedom and the freeze-out temperature. These effects 
are included in our numerical results (for more precise 
analytic treatments see Refs. [DIIIIIIH])- 

We now introduce a simple prescription for comput¬ 
ing the thermal average of the forbidden annihilation 
rate. Detailed balance states that the right-hand side 
of Eq. ® vanishes in equilibrium, = n^. There¬ 
fore, the lorbidden annihilation rate is related to the rate 
of the inverse process, which proceeds at 0 temperature, 
1’) ~ Oidlrn'^d- We find. 
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where A = (m.y^ —m^)/m^ is the relative mass splitting, 
X = m^/T, and /a = (A3/2(2+A)3/2(2+A(2+A)))/(l+ 
A)^. The exponential suppression comes from the form 
of the equilibrium number density for non-relativistic 
species, = g{rnT 12 'kYI'^ exp(—m/T), where 5 = 2 (3) 

for 'i/i ( 7 d), and we have assumed zero chemical potential. 
Note that the approximation of the forbidden cross sec¬ 
tion in Ref. [5] has an incorrect exponential dependence 
on Ax. 

We obtain the forbidden relic density by plugging 
Eq. (§ into Eq. ([^ and integrating the cross section 


from freeze-out to the present in order to account for 
annihilations after freeze-out (see for example Ref. |27)i. 
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where Xf = m^/Tf ~ 10 — 25 and gA{xf) = 

KfA)-\l - is an 0 ( 1 ) 

function. Note that we indicate with the total relic 
density of if) and '0. Eq. (|^ shows that the forbidden relic 
density is exponentially enhanced as A increases. Equiv¬ 
alently, fixing the relic density to the observed value, the 
DM mass is exponentially lighter than the weak scale. 

We show the relic density, as a function of A, in the 
right panel of Fig. Our numerical results here, and 
throughout this letter, utilize Micr0MEGASv4 [53] to solve 
the Boltzmann equations and we have verified that they 
agree with Eq. (|^. The left of the figure, A < 0, corre¬ 
sponds to the conventional case where the relic density is 
too small for light DM masses. As we enter the forbidden 
region, A > 0, the relic density exponentially increases 
until it achieves the correct value. The standard lore is 
that forbidden channels are only relevant in highly de¬ 
generate scenarios, A ^ 1 (this was stated by Ref. | 8 ] 
which implicitly assumes weak scale DM). However, we 
see from Fig. that light DM calls for an 0(1) splitting. 

On the left side of Fig. we show the value of A that 
corresponds to the observed DM abundance, as a func¬ 
tion of the DM mass. For m,/, > 1 MeV, we assume 
that the dark sector is in thermal contact with the SM, 
Tdark = TgM- Lighter masses require DM to be ther¬ 
mally decoupled and cooler, Tdark < Tsm, due to con¬ 
straints on the number of relativistic degrees of freedom 
from Big Bang Nucleosynthesis (BBN) (SO] |3T| and the 
CMB [15]. For < 1 MeV, we adopt a decoupled dark 
sector scenario, consistent with these constraints, that 
we describe below. We find that DM masses down to 
the keV scale are accommodated (DM with a sub-keV 
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Splitting Freeze-Out Temperature 




FIG. 2. The left side shows the splitting A, versus with fixed to the observed abundance. Different curves correspond 
to different values of ad- For > 1 MeV the dark sector temperature equals the SM temperature, while for m-^ < 1 MeV 
the dark sector is cooler than the SM as described in the Thermally Decoupled Dark Sector section. The right side shows the 
freeze-out temperature, xj = as a function of m-^, with A fixed to the value shown to the left. 


mass is too warm, causing problems with structure for¬ 
mation [32H11])- The right side of Fig. shows the value 
of the freeze-out temperature, x/ ^ 10 — 25. 

Self-Interactions: Sizable DM self-interactions may 
leave observable imprints on astrophysical observa¬ 
tions m- The self-interaction cross section is con¬ 
strained by cluster mergers [55MT] and halo shapes m 
1^ . asilm^ < 1 cm^/g ~ 5 x 10 ® MeV A cross 
section near this limit, crsi/rn-tp ^0.1 cm^/g, may help 
resolve the core-cusp HU and too big to fail m prob¬ 
lems of structure formation (see for example Refs. [H- 
m), although baryonic feedback may be important for 
addressing these problems . Self-interactions may 

also be indicated by recent cluster observations [5^ |23] . 
However, an observable self-interacting cross section is 
orders of magnitude larger than the annihilation cross 
section at freeze-out, making self-interactions irrelevant 
for many models of weak scale DM (models with light 
mediators can be an exception [T51148 ) ). 

Forbidden DM naturally achieves a large hierarchy be¬ 
tween the self-interaction cross section, which is unsup¬ 
pressed, and the annihilation rate, which is exponentially 
suppressed after thermal averaging. Self-interactions are 
velocity-independent and are generated by the second di¬ 
agram of Fig. leading to the cross section. 



Signals from Kinetic Mixing: The sizes of direct 
and indirect detection signals depend on how the dark 
sector couples to the SM. We consider the possibility 
that the dark photon kinetically mixes with the SM pho¬ 
ton [T5] . 

( 6 ) 


where any size for e, which characterizes the strength 
of mixing, is technically natural. The kinetic mixing is 
removed (at leading order in e) by shifting the photon, 
—>■ — e 7 ((, inducing a coupling between the dark 

photon and the electromagnetic current, 

The left panel of Fig. shows the plane, fix¬ 

ing ad = 0.1 and choosing the splitting A at each point 
so that the relic density matches the observed value, as 
in Fig. In the lower gray region, the dark sector is 
thermally decoupled from the SM before freeze-out and 
the cosmology depends on the initial conditions. We 
consider this possibility below. Limits on the dark pho¬ 
ton are shown from beam dump experiments |50H54] and 
SN1987A cooling [501115]. 

Kinetic mixing allows DM to annihilate to charged 
states, such as electrons, through off-shell dark photons 
as in the third diagram of Fig. The cross section of 
annihilations into electrons is 
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where /ia = (5 + A(2 -|- A)(5A(2 -l-A) — 

6 ))/((-l +A)2(l +A)4(3 +A)2). For ad « 0.1, 
observable self-interactions are realizable for Forbidden 
DM lighter than ~ 100 MeV, as shown in the left panel 
of Fig. 1^ In the right panel we show the relevant 
parameter space in the {m-,p,ad) plane for a dark sector 
decoupled from the SM. 


This cross section dominates over forbidden annihilations 
in the upper gray region of Fig. By the recombination 
epoch, forbidden annihilations shut off, but annihilations 
into charged states are still active and are constrained 
by CMB observations [IllMISj. The purple region of 
Fig.[^is excluded by Planck m, where we have included 
the efficiency for annihilations to deposit energy into the 
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Signals from Kinetic Mixing (a^=0. 1) Thermally Decoupled Dark Sector (e=0) 



FIG. 3. The left side shows constraints on the kinetic mixing parameter, e, versus m^, with A chosen to match the observed relic 
density and ad = 0.1. Annihilation to SM states dominates over forbidden channels in the upper gray region, and the dark sector 
is thermally decoupled from the SM in the lower gray region. Limits are shown from beam dump experiments (orange) |50H54 | . 
supernovae cooling (blue) |501 155| . Planck (purple) |151156| . and direct detection (pink) | 57H59| . The dashed brown (black) 
curves show the projected reach of SuperCDMS SNOLAB [2] (electron scattering [60]). The red shaded area (dashed curve) 
shows the approximate sensitivity of current observations to DM self-interactions, asijm^ > 1 (0.1) cm^/g |19i 1201 IH5l441| . 
The right side shows constraints on the dark sector when it is thermally decoupled from the SM. In the upper gray area 
A'lp —^ 2-)/) dominates over forbidden annihilations. In the lower gray area DM is overabundant. The red shading and dashed 
curve represent the same values for the self-interaction cross section as in the left panel. In the purple shaded area the DM mass 
is too small to be simultaneously consistent with the Pauli exclusion principle and the densities observed in dwarf spheroidal 
galaxies |35H37| . 


photon plasma from Ref. jSS] . The CMB limit supersedes 
the present reach of diffuse gamma and X-ray observa¬ 
tions [Hi] . 

Kinetic mixing also allows DM to scatter against nu¬ 
clei, as in the fourth diagram of Fig. [l] The DM-nucleon 
cross section is 


cr-ipN 


2 Wnaadn'^p 


( 8 ) 


where = m^mp/{m^ + trip) is the reduced mass of 
DM and the proton. In Fig. we show the strongest 
present limits from direct detection, which, moving from 
heavier to lighter DM mass, come from LUX [^, Super¬ 
CDMS Soudan [^, and CDMSlite [57]. We also show 
the projected sensitivity of SuperCDMS SNOLAB [2], 
which will probe a significant fraction of parameter space. 
DM can also scatter against electrons and we show the 
estimated reach of a future germanium detector [50], al¬ 
though it is superseded in this model by the Planck con¬ 
straint. 

Thermally Decoupled Dark Sector: We now con¬ 
sider the possibility that the dark sector is thermally de¬ 
coupled from the SM during freeze-out, e —>■ 0. Our treat¬ 
ment of the relic density assumes that the dark photons 
remain in equilibrium during freeze-out, with zero chemi¬ 
cal potential, as happens if the dark photons are thermal- 
ized with radiation. In the e —>■ 0 limit, we assume there 
is dark radiation, n, that couples to the hidden photon. 


CD qngdn-f^n-fl^ (9) 


where m„ ^ and ^ 1 is the charge of n under 
the dark force. We assume that (/„ is large enough to 
keep 7 d in equilibrium with n but small enough to pre¬ 
vent ipip —>■ nn decays from dominating over forbidden 
annihilations. A large range of parameters satisfies these 
conditions, 10“^° ^ 9n ^ 10“^^. For < 1 eV, n is a 
warm, subdominant, component of DM that contributes 
less than 10% of the DM energy density, satisfying con¬ 
straints on warm DM [32] . 

In general, the dark sector has a different tempera¬ 
ture than the SM when the two sectors are thermally 
decoupled. We assume that the two sectors begin with 
a common temperature above the weak scale, Tq > v. 
Then, the relative temperatures of the two sectors is de¬ 
termined by the requirement that they separately con¬ 
serve entropy [55], 


Tdark f g^i^jTsM) g^^^HTo) .... 

Tsm \ (To) ^ ^ 

In our model, the hidden sector becomes cooler than the 
SM because more states freeze-out in the SM sector. At 
low temperatures, T m.^, Tjark ~ 0.5Tsm- Because of 
the smaller dark temperature, the hidden sector is con¬ 
sistent with constraints on the number of relativistic de¬ 
grees of freedom from BBN [301 E] and the CMB [15] , 
including when ^ Tbbn ~ 1 MeV. We computed 
these constraints in the presence of a dark Higgs with the 
same mass as the dark photon. 

To the right of Fig. we show the parameter space of 
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the decoupled scenario in the (m^, a^) plane, fixing A to 
match the observed relic density. In the upper gray area, 
> 2m^ and Atp —)■ 2'ijj dominates over forbidden 
annihilations. In the lower gray region, the abundance 
is too large for any choice of A. The shaded red region 
(red dashed line) indicates the limit (approximate observ¬ 
able reach) for self-interactions, crsz/m^ > 1 (0-1) cm^/g. 
Very light fermionic DM, < keV, is not allowed be¬ 
cause the Pauli exclusion principle is inconsistent with 
observed densities of dwarf spheroidal galaxies (purple 
shading) [551 - 157] . A comparable bound is obtained from 
Lyman-a observations [55H55] and also applies to bosonic 
DM. However it can be relaxed if we give up the assump¬ 
tion of a common temperature for the two sectors above 
the weak scale. 

Conclusions: In this letter we have proposed the For¬ 
bidden DM framework. DM may dominantly annihilate 
into heavier states, such that the exponential suppression 
of the thermally averaged cross section allows for DM 
that is exponentially lighter than the weak scale. For¬ 
bidden channels shut off at low temperatures, naturally 
evading CMB constraints. Self-interactions are unsup¬ 
pressed, and potentially observable for Forbidden DM 
lighter than about 100 MeV. We have illustrated For¬ 
bidden DM with a sample model, where fermionic DM 
annihilates into dark photons. We leave consideration of 
more forbidden models for future work, including scenar¬ 
ios where DM annihilates through p-wave processes or is 
a composite state of a strongly coupled hidden sector. 

Acknowledgments: We thank James Bullock, Rou- 
ven Essig, Tracy Slatyer, and Neal Weiner for helpful dis¬ 
cussions. JTR thanks Tel Aviv University for hospitality 
while this work was completed. RTD was supported by 
the NSF grant PHY-0907744. 


* dagnolo@ias.edu 
^ ruderman@nyu.edu 

[1] E. W. Kolb and M. S. Turner, Front. Phys. 69, 1 (1990). 

[2] P. Cushman, C. Galbiati, D. N. McKinsey, H. Robertson, 
T. M. P. Tait, D. Bauer, A. Borgland and B. Cabrera et 
ai, arXiv: 1310.8327 [hep-ex]. 

[3] J. L. Feng and J. Kumar, Phys. Rev. Lett. 101, 231301 
(2008) [arXiv:0803.4196 [hep-ph]]. 

[4] S. Nussinov, Phys. Lett. B 165, 55 (1985). 

[5] D. E. Kaplan, M. A. Luty and K. M. Zurek, Phys. Rev. 
D 79, 115016 (2009) [arXiv:0901.4117 [hep-ph]]. 

[6] E. D. Carlson, M. E. Machacek and L. J. Hall, Astrophys. 
J. 398, 43 (1992). 

[7] Y. Hochberg, E. Kuflik, T. Volansky and J. G. Wacker, 
Phys. Rev. Lett. 113, 171301 (2014) [arXiv: 1402.5143 
[hep-ph]]. 

[8] K. Griest and D. Seckel, Phys. Rev. D 43, 3191 (1991). 

[9] G. B. Jackson, G. Servant, G. Shaughnessy, T. M. P. Tait 
and M. Taoso, JCAP 1004, 004 (2010) [arXiv:0912.0004 
[hep-ph]]. 

[10] S. Tulin, H. B. Yu and K. M. Zurek, Phys. Rev. D 87, 


no. 3, 036011 (2013) [arXiv: 1208.0009 [hep-ph]]. 

[11] G. B. Jackson, G. Servant, G. Shaughnessy, T. M. P. Tait 
and M. Taoso, JCAP 1307, 021 (2013) [arXiv: 1302.1802 
[hep-ph]]. 

[12] C. B. Jackson, G. Servant, G. Shaughnessy, T. M. P. Tait 
and M. Taoso, JCAP 1307, 006 (2013) [arXiv: 1303.4717 
[hep-ph]]. 

[13] M. Pospelov and J. Pradler, Phys. Rev. D 82, 103514 
(2010) [arXiv:1006.4172 [hep-ph]]. 

[14] N. Padmanabhan and D. P. Finkbeiner, Phys. Rev. D 
72, 023508 (2005) [astro-ph/0503486]. 

[15] P. A. R. Ade et al. [Planck Collaboration], 
arXiv: 1502.01589 [astro-ph.CO]. 

[16] D. N. Spergel and P. J. Steinhardt, Phys. Rev. Lett. 84, 
3760 (2000) [astro-ph/9909386]. 

[17] W. J. G. de Blok, Adv. Astron. 2010, 789293 (2010) 
[arXiv:0910.3538 [astro-ph.CO]]. 

[18] M. Boylan-Kolchin, J. S. Bullock and M. Kapling- 
hat, Mon. Not. Roy. Astron. Soc. 415, L40 (2011) 
[arXiv: 1103.0007 [astro-ph.CO]]. 

[19] M. Rocha, A. H. G. Peter, J. S. Bullock, M. Kapling- 
hat, S. Garrison-Kimmel, J. Onorbe and L. A. Mous- 
takas, Mon. Not. Roy. Astron. Soc. 430, 81 (2013) 
[arXiv: 1208.3025 [astro-ph.CO]]. 

[20] A. H. G. Peter, M. Rocha, J. S. Bullock and M. Kapling- 
hat, Mon. Not. Roy. Astron. Soc. 430, 105 (2013) 
[arXiv: 1208.3026 [astro-ph.CO]]. 

[21] J. Zavala, M. Vogelsberger and M. G. Walker, Monthly 
Notices of the Royal Astronomical Society: Letters 431, 
L20 (2013) [arXiv:1211.6426 [astro-ph.CO]]. 

[22] R. Massey, L. Williams, R. Smit, M. Swinbank, 
T. D. Kitching, D. Harvey, M. Jauzac and H. Israel 
et ai, Mon. Not. Roy. Astron. Soc. 449, 3393 (2015) 
[arXiv: 1504.03388 [astro-ph.CO]]. 

[23] F. Kahlhoefer, K. Schmidt-Hoberg, J. Kummer and 
S. Sarkar, arXiv:1504.06576 [astro-ph.CO]. 

[24] M. Pospelov, A. Ritz and M. B. Voloshin, Phys. Lett. B 
662, 53 (2008) [arXiv:0711.4866 [hep-ph]]. 

[25] E. Izaguirre, G. Krnjaic, P. Schuster and N. Toro, 
arXiv:1505.00011 [hep-ph]. 

[26] R. T. D’Agnolo and J. T. Ruderman, to appear. 

[27] P. Gondolo and G. Gelmini, Nucl. Phys. B 360, 145 
(1991). 

[28] G. Steigman, B. Dasgupta and J. F. Beacom, Phys. Rev. 
D 86, 023506 (2012) [arXiv: 1204.3622 [hep-ph]]. 

[29] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, 
arXiv: 1407.6129 [hep-ph]. 

[30] R. H. Gyburt, B. D. Fields, K. A. Olive and E. Skillman, 
Astropart. Phys. 23, 313 (2005) [astro-ph/0408033]. 

[31] R. Gooke, M. Pettini, R. A. Jorgenson, M. T. Murphy 
and C. G. Steidel, arXiv: 1308.3240 [astro-ph.CO]. 

[32] M. Viel, J. Lesgourgues, M. G. Haehnelt, S. Matarrese 
and A. Riotto, Phys. Rev. D 71, 063534 (2005) [astro- 
ph/0501562]. 

[33] A. Loeb and M. Zaldarriaga, Phys. Rev. D 71, 103520 
(2005) [astro-ph/0504112]. 

[34] D. Hooper, M. Kaplinghat, L. E. Strigari and 

K. M. Zurek, Phys. Rev. D 76, 103515 (2007) 

[arXiv:0704.2558 [astro-ph]]. 

[35] S. Tremaine and J. E. Gunn, Phys. Rev. Lett. 42, 407 
(1979). 

[36] A. Boyarsky, O. Ruchayskiy and D. lakubovskyi, JGAP 
0903, 005 (2009) [arXiv:0808.3902 [hep-ph]]. 

[37] S. Horiuchi, P. J. Humphrey, J. Onorbe, K. N. Abazajian, 



6 


M. Kaplinghat and S. Garrison-Kimmel, Phys. Rev. D 
89, no. 2, 025017 (2014) [arXiv:1311.0282 [astro-ph.CO]]. 

[38] M. Markevitch, A. H. Gonzalez, D. Clowe, A. Vikhlinin, 
L. David, W. Forman, C. Jones and S. Murray et al, 
Astrophys. J. 606, 819 (2004) [astro-ph/0309303]. 

[39] D. Clowe, A. Gonzalez and M. Markevitch, Astrophys. 
J. 604, 596 (2004) [astro-ph/0312273]. 

[40] S. W. Randall, M. Markevitch, D. Clowe, A. H. Gon¬ 
zalez and M. Bradac, Astrophys. J. 679, 1173 (2008) 
[arXiv:0704.0261 [astro-ph]]. 

[41] D. Harvey, R. Massey, T. Kitching, A. Taylor and E. Titt- 
ley, Science 347, no. 6229, 1462 (2015) [arXiv: 1503.07675 
[astro-ph.CO]]. 

[42] F. Governato, C. Brook, L. Mayer, A. Brooks, G. Rhee, 
J. Wadsley, P. Jonsson and B. Willman et al, Nature 
463, 203 (2010) [arXiv:0911.2237 [astro-ph.CO]]. 

[43] A. Pontzen and F. Governato, Mon. Not. Roy. Astron. 
Soc. 421, 3464 (2012) [arXiv:1106.0499 ]astro-ph.CO]]. 

[44] F. Governato, A. Zolotov, A. Pontzen, C. Christensen, 
S. H. Oh, A. M. Brooks, T. Quinn and S. Shen et 
al, Mon. Not. Roy. Astron. Soc. 422, 1231 (2012) 
[arXiv:1202.0554 ]astro-ph.CO]]. 

[45] A. Zolotov, A. M. Brooks, B. Willman, F. Governato, 
A. Pontzen, C. Christensen, A. Dekel and T. Quinn et 
al, Astrophys. J. 761, 71 (2012) [arXiv: 1207.0007 [astro- 
ph.CO]]. 

[46] J. L. Feng, M. Kaplinghat and H. B. Yu, Phys. Rev. Lett. 
104, 151301 (2010) ]arXiv:0911.0422 [hep-ph]]. 

[47] M. R. Buckley and P. J. Fox, Phys. Rev. D 81, 083522 
(2010) [arXiv:0911.3898 [hep-ph]]. 

[48] A. Loeb and N. Weiner, Phys. Rev. Lett. 106, 171302 
(2011) [arXiv:1011.6374 [astro-ph.CO]]. 

[49] B. Holdom, Phys. Lett. B 166, 196 (1986). 

[50] J. D. Bjorken, R. Essig, P. Schuster and N. Toro, Phys. 
Rev. D 80, 075018 (2009) ]arXiv:0906.0580 [hep-ph]]. 


[51] R. Essig, R. Harnik, J. Kaplan and N. Toro, Phys. Rev. 
D 82, 113008 (2010) [arXiv: 1008.0636 [hep-ph]]. 

[52] J. Blumlein and J. Brunner, Phys. Lett. B 701, 155 
(2011) ]arXiv:1104.2747 [hep-ex]]. 

[53] S. N. Gninenko, Phys. Lett. B 713, 244 (2012) 
[arXiv: 1204.3583 [hep-ph]]. 

[54] S. Andreas, C. Niebuhr and A. Ringwald, Phys. Rev. D 
86, 095019 (2012) [arXiv: 1209.6083 [hep-ph][. 

[55] J. B. Dent, F. Ferrer and L. M. Krauss, arXiv: 1201.2683 
[astro-ph.CO]. 

[56] T. R. Slatyer, arXiv:1506.03811 [hep-ph]. 

[57] R. Agnese et al. [SuperCDMS Collaboration], Phys. 
Rev. Lett. 112, no. 4, 041302 (2014) ]arXiv: 1309.3259 
[physics.ins-det]]. 

[58] D. S. Akerib et al. [LUX Collaboration], Phys. Rev. Lett. 
112, 091303 (2014) [arXiv:1310.8214 [astro-ph.CO]]. 

[59] R. Agnese et al. [SuperCDMS Collaboration], Phys. Rev. 
Lett. 112, no. 24, 241302 (2014) [arXiv:1402.7137 [hep- 
ex]]. 

[60] R. Essig, J. Mardon and T. Volansky, Phys. Rev. D 85, 
076007 (2012) ]arXiv:1108.5383 [hep-ph]]. 

[61] T. R. Slatyer, N. Padmanabhan and D. P. Finkbeiner, 
Phys. Rev. D 80, 043526 (2009) ]arXiv:0906.1197 [astro- 
ph.CO]]. 

[62] S. Galli, T. R. Slatyer, M. Valdes and F. locco, Phys. Rev. 
D 88, 063502 (2013) [arXiv: 1306.0563 [astro-ph.CO]]. 

[63] M. S. Madhavacheril, N. Sehgal and T. R. Slatyer, 
Phys. Rev. D 89, 103508 (2014) ]arXiv:1310.3815 [astro- 
ph.CO]]. 

[64] R. Essig, E. Kuflik, S. D. McDermott, T. Volansky and 
K. M. Zurek, JHEP 1311, 193 (2013) ]arXiv: 1309.4091 
[hep-ph]]. 

[65] J. L. Feng, H. Tu and H. B. Yu, JCAP 0810, 043 (2008) 
[arXiv:0808.2318 ]hep-ph]]. 



